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CLAIMS 

We claim: 

1 . A method for enhancing an image containing endmembers, comprising: 
selecting at least one measure of disparity; 
5 partitioning at least one first set of said endmembers in said image into at least 

one second set of said endmembers, each said second set incorporating at least 
one site having some spectral content, 
wherein said at least one site is a generic element of a lattice, and 

wherein spatial consistency is imposed on said spectral content of each said site so that 
10 each said at least one second set is associated with a spatially smooth region in said 

image; and ' 

using said at least one measure of disparity, applying a linear mixing model to 
said at least one second set to globally label said at least one second set, 
wherein said method yields improvements in: 
15 assessment of the type and amount of materials in a scene; 

unsupervised clustering of a scene, and 

post-processing smoothing operations suitable for application to diverse 
techniques for terrain categorization and classification. 

20 2. The method of claim 1 in which a multi-grid Gibbs-based algorithm is used to partition 
said at least one first set of said endmembers, 

wherein said algorithm is used to estimate an underlying and imobserved process, X**, a 
discrete labeling Markov Random Field (MRF) process, that associates a label with each 
said site, and 

25 wherein said is used to condition a spectral mixing process (SMP), and 

wherein said X^ identifies each said site that may be treated as homogeneous, and 
wherein said partitioning may proceed with multiple material types incorporated within 
each said site. 

30 3. The method of claim 2 in which said partitioning using said niulti-grid Gibbs-based 
algorithm proceeds in stages from coarse toward fine grid resolutions. 
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wherein the number of said stages is defined by the lattice structure of a neighborhood, 
forming a simple pattern such that a neighborhood of each said lattice element 

comprises near, intermediate, and far neighbors that are specific multiples of the grid 
resolution, a, and 

5 wherein the near neighbors are closest said lattice elements that may touch said site even 

if only at said site's comers, and 

wherein the next 16 said lattice elements are intermediate neighbors, and 
wherein the far neighbors are the 24 said lattice elements on the outside of said next 16 
lattice elements, thus forming perimeter layers away firom said site of interest at multiples 
10 of said G, and 

wherein use of said intermediate and far neighbors facilitates a faster global solution than 
using just said near neighbors, and 

wherein use of said intermediate and far neighbors also enables said algorithm to 
remember label assignments from prior said stages, and 
15 wherein said a determines the spatial sampling of said algorithm at some specific said 

stage of multi-grid processing, such that at full resolution, i.e., a = 1, said site 
corresponds to a single said lattice element, and 

wherein all coarse resolutions involve neighborhoods of said lattice elements forming a 
square in which said a is the number of pixels along one side of said square. 

20 

4. The method of claim 2 in which said Gibbs-based algorithm is initialized conventionally, 
wherein said initialization enables a post-processing smoothing algorithm that may 
operateonthe outputof at least one classification technique. ' 

25 5. The method of claim 2 in which said partitioning is initialized randomly, 

wherein said randomly initialized partitioning algorithm is best implemented as a 
multigrid process using an extended neighborhood system to attain reasonable global 
labeling accuracy,. 
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6. The method of claim 2 in which said partitioning is initiaUzed using a pre-processing 
scheme, 

7. The method of claim 6 in which said pre-processing scheme is a recognized classification 
5 technique. 

8. The niethod of claim 7 in which said recognized classification technique is supervised. 

9. The method of claim 7 in which said recognized classification technique is unsupervised. 

10 

10. The method of claim 2 in which said multi-grid Gibbs-based algorithm uses an energy 
function to model distinctions in said image, 

wherein, using the Gibbs Equivalent Theorem with an appropriately defined graph, 
{S^^\^^^'^^} , said Gibbs-based algorithm may be derived, and 

15 wherein a maximum a posteriori (MAP) estimate may be computed that maximizes 

Vr(J^ I g) by iteratively sampling from the local Gibbs distribution pertaining to each said 
site, s eS^/\ and 

wherein computational intensity is reduced and said global labeling is improved. 

20 11. The method of claim 10 in which said MAP estimate of Jf^, denoted asX/*, is 
calculated ill a first step, termed the Expectation Step, by: 

X''=argmax Pr(X''\X^) and 
the maximum likelihood estimate (MLE) of a vector of proportions on a label 
lattice, p, deiioted as , is calculated in a second step, termed the Maximization Step, by: 

25 )ff = argmaxZ()ff|X\Xi') 

wherein said MAP estimate of X^is used as if it were observed, and 
wherein complex integration over all possible realizable partitions is avoided by said two- 
step process. 
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12. The method of claim 10 in which said MAP estimate is computed using ratios of the 
relationship: Pr(JC^ =x^\x^. = x^,r^s)^FT^^ 



where: ^ and 



Us{x^ ,g) = energy interaction of site s e S^/^ with the neighborhood 



13. The method of claim 1 in which said at least one measure of disparity is selected from the 
group consisting essentially of: spectral angle, Euclidian distance, classical Kolmogorov- 
Smimdv measiu-es, mean adjusted Kolmogorov-Smimov measures, and conibinations 

10 thereof, 

wherein said image is defined with respect to a local neighborhood. 

14. The method of claim 13 in which said local neighborhood comprises at least one non- 
traditional and locally-extended neighborhood. 

15 ' ■ 

15. The method of claim 1 in which parameters of said linear mixing model are estimated 
based on said endmembers in said at least one second set being associated with 
appropriate said sites therein to enable computation of an enhanced spectral mixing 
process (SMP), 

20 wherein the set of material spectra, H, is represented as 

and ciH is the set of material spectra at a pixel site s, Ngnds is the total number of 

said endmembers, and the spectra for the k^^ said endmember is given by 
/(W = G A) , if A^*^ G /f, , 



16, The method of claim 1 further employing a simulated annealing method of partitioning 
hyperspectral imagery, said simulated annealing inethod initialized by a supervised 
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classification method to provide spatially smooth class labeling for terrain mapping 
applications. 

17. The method of claim 1 providing an estimate of said model as a Gibbs distribution 
^ defined over a symmetric spatial neighborhood system that is based on an energy 

fiinction characterizing spectral disparities in both Euclidean distance and spectral angle, 
wherein said energy fiinction depends on distance measures termed disparity metrics that 
provide a measure of the dissimilarity between an individual said site, s, and said 
individual site's neighbor, r. 

18. An efficient and accurate method for extracting features fi'om hyperspectral data 
representing a scene by implementing a spectrally-optimal supervised classification 
algorithm that has been smoothed by a post-processing routine that imposes 
spectral/spatial constraints defined by the Gibbs prior probability distribution, Pr(^). 

19. A method for extracting features fi*om hyperspectral data representing a scene by 
initializing a partitioning algorithm with the results of a classifier to improve 



classification by providing initial estimates based on labels that, as applied under the 



proper conditions, are spectrally optimal, comprising: 

a) performing a classification of the scene using a supervised classifier that selects 

the labeling Xs at each said site s that maximizes the probability Pr (^j' ) 

of a process that is approximately of the form of 



where: 




= the mean of the spectral vectors for the class label, xf 




= the covariance of said spectral vectors; and 



b) initializing by the results in Step a), performing a spectral/spatial 
partitioning of said scene using the relationship: 



PTix:=x:\x:=x:,r^s)=Mx:=x:\x^ =x:,r^^r')-Y'^"''^''' 



s 
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where: , 



25 



Us{ ^) = energy interaction of site s e S^/^ with the neighborhood ^ . 

5 20. A method employing a simulated anneaHng method of partitioning hyperspectral imagery 
for extracting features from hyperspectral data representing a scene, comprising: 

using a Bayesian framework to develop a 2-step supervised Gibbs-based 
classification algorithm, 
wherein said algorithm is capable of performing high quality spatially smooth labeling of 
10 said hyperspectral data; and 

using a linear classifier to initialize said Gibbs-based partitioning algorithm, 
wherein said initializing results in improved label accuracy and smoothness as compared 
to using only said linear classifier without said Gibbs-based partitioning algorithm, and 
wherein global labeling accuracy is increased as compared to a stand-alone randomly 
15 initialized said Gibbs-based partitioning algorithni, and 

wherein said initializing said' Gibbs-baseid partitioning algorithni with said linear 
classifier also reduces the computation by eliminating the need for a mlilti-grid process 
and by allowing said simulated annealing to start at a cooler temperature. 
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